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Abstract 

The study of the three-body problem with short-range attractive two-body forces has a rich 
history going back to the 1930's. Recent applications of effective field theory methods to atomic and 
nuclear physics have produced a much improved understanding of this problem, and we elucidate 
some of the issues using renormalization group ideas applied to precise nonperturbative calculations. 
These calculations provide 11-12 digits of precision for the binding energies in the infinite cutoff 
limit. The method starts with this limit as an approximation to an effective theory and allows 
cutoff dependence to be systematically computed as an expansion in powers of inverse cutoffs 
and logarithms of the cutoff. Renormalization of three-body bound states requires a short range 
three-body interaction, with a coupling that is governed by a precisely mapped limit cycle of the 
renormalization group. Additional three-body irrelevant interactions must be determined to control 
subleading dependence on the cutoff and this control is essential for an effective field theory since 
the continuum limit is not likely to match physical systems {e.g., few-nucleon bound and scattering 
states at low energy). Leading order calculations precise to 11-12 digits allow clear identification 
of subleading corrections, but these corrections have not been computed. 

PACS numbers: ll.lO.Gh, 03.65..Ge, 05.10.Cc, 21.45. +v, 36.40.-c 
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I. INTRODUCTION 



It has been known since the 1930's that a non-relativistic three-body system with short- 
range two-body potentials has a pecuhar behavior, with a binding energy that is unex- 
pectedly large. In fact the binding energy becomes infinite if the range of the attractive 
two-body potential is shrunk to zero Q. For over 30 years it has also been known that the 
three-body bound state spectrum exhibits an Efimov effect 05l3- There are an increasing 
number of three-body bound states as the two-body effective range is reduced or as a high 
momentum cutoff is increased, as long as the two-body bound state energy is held fixed. 
For sufficiently large cutoff. A, the ratio of the energies of two successive three-body bound 
states approaches an analytically fixed limit of approximately 515.035 for energies much 
deeper than the two-body binding energy j^]. 

More recently, Bedaque, Hammer, and van Kolck have shown that the short-range three- 
body problem is renormalizable, but only if a point-like three-body interaction is added 
0,0, 01. The dimensionless version of the three-body coupling strength [which they denote 
by H{A) but we denote by G3(A)] has an unexpected dependence on A, namely G3(A) 
has a periodic dependence on In A. In each period the value of G3(A) steadily increases 
until it reaches plus infinity. It jumps to minus infinity and then steadily increases again. 
This peculiar behavior provides a rare example of a renormalization group limit cycle |^ 
lol |. with the exciting possibility that scaling behavior near a limit cycle might be studied 



experimentally [11]. 

Unfortunately, the two experimental candidates to be Efimov-like, namely the triton and 
the atomic helium trimer, do not exhibit the infinite set of shallow Efimov states because 
that would require the nucleon-nucleon S-waves or helium dimer systems to have a bound 
state precisely at threshold |2,13]- This is not the case experimentally. The reality is that the 
triton has only one bound state. In the case of the helium trimer numerical calculations for 
realistic potential models indicate that the trimer has two bound states. These departures 
from the infinite Efimov limit raises the question: what are the corrections to the Efimov 
limit in the three-body system when the two-body system does not have a bound state 
precisely at threshold? How large are the corrections in the case that there are only one or 
two Efimov states rather than an infinite number? 

Ideally, one would answer this question for the most plausible three-body Hamiltonians 
that have been proposed for describing the triton or the helium trimer [see refs. in But 
realistic three-body Hamiltonians are very demanding to study. Their three-body eigenfunc- 
tions are functions of three variables in the simplest case (S-states), making both analytical 
and numerical analysis exceptionally difficult (although far from impossible). Fortunately, 
it is possible to formulate a cutoff form of a two-body zero range potential (a delta-function 
potential V"(r) = g6^{r)) in which the Schrodinger equation wave function can be expressed 
exactly in terms of a reduced wave function 0(r) that depends only on one radial variable 
rather than three variables. The Schrodinger equation itself becomes a one-dimensional in- 
tegral equation for 0(r) and the energy eigenvalue E (although with auxiliary functions that 
involve integrals over known functions). 

The one- dimensional integral equation can be solved numerically and quickly to machine 
precision (around 12 decimal places in our calculations) on a PC, although achieving this 
high precision involves careful attention to details of how the reduced wave function and 
the integral equation are discretized. High numerical accuracy is invaluable for testing the 
validity of analytic but approximate formulae for the corrections to the Efimov limit. 
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The integral equation will be formulated and solved for the Fourier transform of 
the reduced wave function. Uniformly valid expansions {e.g., with uniform convergence to 
the exact solution over the entire range of momenta) are employed to systematically dissect 
the bound state equation so that high-momentum behavior is isolated analytically. Low, 
intermediate and high momentum scales in the bound-state wavefunction are isolated by 
defining approximations to the full wavefunction that are valid in each of these three regions, 
and then showing that these can be assembled to produce an approximation that is valid for 
all momenta. The wave-functions valid at low and high momenta "communicate" through 
a single phase found in an analytic solution for the intermediate momenta approximation, 
all to leading order. 

The outcome of the analytic and numerical analysis of this article is that Efimov three- 
body bound states have energy eigenvalues whose infinite cutoff limits can be precisely 
computed, and a simple renormalization group analysis suggests the leading corrections are 
of order a/B2 ln(A)/A and B^^^ ln(A)/A^ relative to -83"^ itself. B2 is the binding energy 
of the two-body state near threshold. But these leading corrections have coefficients which 
depend on the details of the two- and three-body interactions in the Hamiltonian, such as 
the detailed shape of the two-body cut off potential, which suggests that the actual size of 
corrections to the Efimov limit for the triton or helium trimer could depend on the details 
of the Hamiltonians that best characterize these systems. 

There is a second reason for writing this article. The Hamiltonian studied in this arti- 
cle exhibits renormalization group limit cycle behavior, at least to the numerical accuracy 
achieved. We develop a renormalization group description of the Hamiltonian that takes 
into account the leading corrections to the Efimov limit as derived from the integral equa- 
tion. The renormalization group description takes the form of an extended Gell-Mann-Low 
analysis jl^ involving two-body and three-body coupling parameters. The extended Gell- 
Mann-Low analysis provided here accounts for the limit cycle and for the leading corrections 
that are of relative order a/B^ ln(A)/A as well as corrections of order \fB^j K. This requires 
two two-body couplings and one three-body coupling. The analysis reported here is some- 
what similar to the analyses produced in effective field theory (EFT). However the analysis 
offered here differs from EFT in both the two- and three-body sectors because it starts from 
a non-perturbative fixed point in the two-body sector and from the three-body limit cycle 
in the three-body sector. There are no new outcomes in our approach that differ from the 
EFT analysis for the two-body sector. It is only in the three-body sector that, due to the 
limit cycle, the renormalization group analysis provides a more straightforward approach 
than anything published to date using EFT. 

A major concept of the renormalization group is the concept of universality. The concept 
of universality was developed for the case of renormalization group fixed points and has 
two parts. The first part (due mostly to Franz Wegner jl^) is an analysis of infinitesimal 
departures from the fixed point Hamiltonian, such as infinitesimal departures from the fixed 
point of the two-body Hamiltonian. Franz Wegner also formulated a classification of op- 
erators that are relevant, marginal or irrelevant. This classification requires adjustments 
before it applies to infinitesimal departures from a renormalization group limit cycle. One 
adjustment is that whenever there is a limit cycle, there is an infinitesimal marginal operator 
whose role is to correspond to an infinitesimal shift on the limit cycle itself. The coefficient 
of this "limit cycle shift" operator is an angular variable 6', and all physics is invariant to 
changes in Q by I-r . Other adjustments are more minor. 

We will present our analysis in a number of stages, with analysis of the integral equation 
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postponed until the last two stages. The stages are as follows: 

• The two-body renormalization group, neglecting irrelevant operators (Section II). 

• The two-body renormalization group with the leading irrelevant operator (Section III). 

• The three-body renormalization group, neglecting irrelevant operators, and the possi- 
bility of a limit cycle (Section IV). 

• The three-body renormalization group with the leading irrelevant operator (Section V). 

• Derivation of the integral equation in the presence of a Gaussian-like cutoff (Sec- 
tion VI). 

• Methods of solution of the integral equation: method for analysis of limiting behavior 
for large A, applied to a simplified example (Section VII). 

• A renormalized equation in the three-body case (Section VIII). 

• Discretization of the integral equations with exponentially small errors (Sections IX 
and X). 

• Analytic and numerical results (Section XI). 

The initial discussions of the three-body renormalization group equations draw on results 
reported later, so a full understanding of this work probably requires two passes. Additional 
details and precise calculations can be found in Richard Mohr's thesis jl^. For a detailed 
account of previous work on the three-body system with short-range interactions, see Ref. 



Ill and references therein. 



II. THE TWO-BODY RENORMALIZATION GROUP 

We develop the renormalization group equations first for the two-body sector. These 
equations serve unchanged as a subset of the renormalization group equations for the three- 
body sector. In its simplest form, excluding irrelevant operators, the two-body renormal- 
ization group equation is for a single dimensionless coupling constant 6*2 (A). The equation 
for G2(A) has a fixed point solution at a value G2 for which the two-body binding energy is 
zero. The fixed point solution for G2(A) provides a contrast to the limit cycle solution for 
6*3 (A), and it is interesting in its own right because it is a non-free fixed point far removed 
from the free fixed point [l5|. The existence of this fixed point leads to scaling behavior that 
seems unnatural if one seeks to expand around the free theory Its existence brings 



into question the use of free operators (i.e., powers of fields and derivatives with dimension 
given by a free scaling analysis) instead of operators with good scaling behavior (i.e., eigen- 
operators of a linearized RG transformation about the interacting fixed point 13].) For a 
complete RG analysis using a slightly different formalism, see (l5| . 

The bound state equation and the scattering amplitude for a separable two-body Hamil- 
tonian are known; we report the result as the starting point for analysis. We consider a 
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two-body Hamiltonian with two identical Bose particles of mass M = 1/2 in the center-of- 
mass frame. This mass is chosen to simplify the algebra in the three-body problem. The 
equivalent one-body Hamiltonian is: 

H^2p^ + V , (1) 
where we choose the separable two-body potential V: 

y(p,p') = -^c/A(p)c/A(p'). (2) 

Ua{p) can in principle be any real function of the three- momentum p and the cutoff A that 
goes to zero as p becomes much larger than A and to one as p goes to zero, but symmetries 
should be respected if possible. We have written the interaction using a dimensionless 
coupling, G2, because we need to disentangle the scaling dependence of couplings that 
serve as coordinates in a space of operators from the scaling dependence of the operators 
themselves. It is much simpler to develop RG-improved perturbation theory about a non- 
free fixed point using the parameters that must remain small near the fixed point, and these 
are the dimensionless couplings. 

The exact bound state equation reduces to 

where the integration volume is dr = d^p/ (27r)^. The K- matrix is equivalent to the T-matrix 
but uses standing wave boundary conditions that produce a principal value prescription for 
scattering integrals. It satisfies the integral equation 

.,P,p';.) = nP,P').P/|^M^. 

For a separable potential this equation is easily solved: 

;.-(p,p.^)._iM^, (5) 

where 

Once again, this is an exact result, which allows us to identify and explore the region about 
the fixed point. 

The K-matrix is directly related to the S-wave phase shift and thereby to the effective 
range expansion for low energy scattering using 

Pcotft) = -— = + + (7) 

where K[E) is the on-shell K-matrix with p = p',E = 2p^, a is the scattering length and is 
the effective range; and we will confine our attention to S-waves. This relationship is ideally 
suited to the development of a Gell-Mann-Low analysis. We take advantage of the fact that 
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we can solve this two-body problem exactly in the limit A ^ oo, using a Gell-Mann-Low 
analysis to identify all operators required, rather than attempting an exact Wilsonian RG 
analysis 

This limit yields (see for example jiol ]) 

p cot{5p) = , D{E) = 0. (8) 

The scattering length is infinite, so there is a bound state at zero energy. A binding energy 
can only be introduced by adding a scale, which separates infrared and ultraviolet scaling 
behavior and which therefore must be associated with moving away from an ultraviolet fixed 
point. We do not provide a full analysis, but for any couplings that have a continuum hmit, 

p cot(5p) = -^/B^. (9) 

This is true for any momentum and any cutoff because an exact RG produces 'physical' 
results that are completely independent of the cutoff and which therefore must be the same as 
these results in the limit where the cutoff is taken to infinity. This can be used to understand 
the limits of a Gell-Mann-Low analysis. In general, one adjusts the irrelevant operators in 
an EFT Hamiltonian so that the effective range expansion matches data. This requires 
these couplings to be tuned away from their continuum RG trajectory, and we know that 
running the transformation backward will cause irrelevant deviations from continuum theory 
to deviate exponentially from continuum trajectories in the ultraviolet limit. At some point 
the irrelevant couplings start to become large and the RG-improved perturbative analysis 
breaks down; the theory becomes unnatural in the parlance of effective field theory. There 
is also an issue of how rapidly an expansion of the actual irrelevant operators in powers of 
free operators converges, but we do not investigate this issue here and turn instead to the 
use of a smooth Gaussian cutoff. 

The Hamiltonians we want to study must be approximated as deviations from the fixed 
point, and in general we will approximate the fixed point and the relevant and irrelevant 
operators using 

U^ip) = (1 + /i2(A) p^ + ...) exp(-pVA')- (10) 

We can adjust 6*2 (A) to fix the binding energy by solving Eq. 0, but we see that it depends 
on h2 can be fixed by also insisting that the effective range be zero, which is true for any 
continuum theory and introduces the first place at which effective theories must be allowed 
to deviate from the continuum limit for practical application. This requires the introduction 
of irrelevant couplings that become of order 1 when A = l/rg if the effective range, r^, is 
not zero, as we will see below. 

It is clearly quite easy to complete a non-perturbative analysis of the two-body problem, 
but we will use the Gaussian cutoff and analyze the problem using the method of uniformly 
valid expansions, which is a key tool in our three-body analysis that is presented in Section 
VII. We do not outline the method here but only point out its appearance, saving the 
discussion of details. We will first approximate the fixed point and relevant operator for the 
fixed point using 

Ua{p) = exp(-pVA2) . (11) 

The equation we seek is an equation for dG2/dA, expressed as a function of G2 and A. 
Since G2 is dimensionless, while A is not, the structure of the equation must be 

rlG 

^-7^ = Pg,{G2) . (12) 
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where (io^iGT) is a function to be determined, in this case from the equation for 1/^2- 
We assume there is a fixed point Gg, and the function ^02(^2) can be determined as an 
expansion in powers of G2 — G^2- The result takes the form: 



with 



Gl = 4(27r)=^/^ and 



C2 = 



- GlY + 
(2-vr) 



(13) 



(14) 



2(27r)5/2 ■ 

Note that there is no exact fixed point solution when a Gaussian cutoff is used, so the 
value of couplings at the fixed point will change in our Gell-Mann-Low analysis as we add 
couplings. We will provide only the beginning of the analysis that leads to RG equations. 
The equation for I/G2 can be rewritten as 



A 

G'2 



27r2 



dp 



p2exp(-2pVA2) 
{2p^ + B2) 



(15) 



where we can employ a uniformly valid expansion at next-to-leading order to write 



A _ 1 /• 



dp 



p 



2p2 



1 



2p2 + B2 



A^ +2^"^^ 



-V/A^) 1 - S 




~2p^) 

(16) 

The first term is the leading order approximation to the integrand in the region of integration 
where p ^ B2 <^ A, the second term is the approximation in the region where -B2 <^ P ~ A, 
and the third term is the approximation in the intermediate region which cancels the leading 
errors in the other two approximations in the regions where they are not valid. This division 
is not necessary here, where we know the integrand, but it immediately leads to: 



1 

G'2 



1 

G* 



+ 



Bo 



Sv^TT A 2(27r)3/2A 



+ 




(17) 



This equation makes it clear that our potential can be used to approximate both the fixed 
point and the relevant operator for that fixed point, where our approximate relevant operator 
is used to control the binding energy {i.e., the scattering length). 



III. TWO-BODY RENORMALIZATION GROUP WITH AN IRRELEVANT OP- 
ERATOR 

We now expand the renormalization group equations in the two-body sector to include 
a coupling coefficient for the leading irrelevant operator. The purpose of this exercise is to 
identify and control the leading cutoff-dependent correction in the effective range expansion. 
We introduce a new coupling into the two-body Hamiltonian that, when varied, enables this 
correction to be held fixed even as the cutoff A is replaced by an effective cutoff A'. 

We want to contrast the form of the leading irrelevant operator in the two-body sector 
with the leading irrelevant operator obtained in perturbation theory. Thus, as a preliminary, 
we look at irrelevant operators in first-order perturbation theory for small G2- 

To first order, the K-matrix is 

i^(p,p') = ^C/A(p)C/A(p')- (18) 
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The dominant terms in this expression, for fixed momenta and large A, are 



;.-(p,pO.^M^-M±p!). (19) 

In this case, the leading irrelevant terms are the terms involving and p'^; to keep K 
constant while changing A, one could generalize the function ^/a to become 

[/a(p) =[l + /^2(A)^) exp(-pVA2) . (20) 

Then the renormalization group equations needed to hold K fixed to first order in an ex- 
pansion in powers of p^/A^ are: 

and 

^IK - ^2(A)]) = 3G2(A)[1 - /12(A)] . (22) 

Now we repeat this analysis for the full nonperturbative amplitude. We identify the 
leading nonperturbative A-dependent correction to the scattering amplitude near the fixed 
point in G2- It is sufficient to identify the leading correction to the function D{E). We shall 
assume that E is of order B2. We can rewrite the equation for D{E), Eq. (6), using Eq. (jH} 
to give: 

uUp) K(v) 



^'^)-^/. *|(VTly-(#^^l ■ '''' 

This integral is finite as A ^ 00, so it should have an expansion in inverse powers of A, 
although we must allow for the possibility that logarithms of the cutoff will appear. We 
can make some progress before specifying Uf^ by using a uniformly valid expansion for the 
integrand (as discussed in Section VII) and the constraint that the cutoff function goes to 
zero for large momenta and to one for small momenta. To leading order the uniformly valid 
expansion gives us: 

2^- ^Jo ^^i(2p2 + i?2)(2p2-E)+^" Vl • ^^^^ 

Again, the first term is the leading order approximation to the integrand in the region of 
integration where p B2, E <^ A, the second term is the approximation in the region where 
B2, E p ^ A, and the third term subtracts the approximation in the intermediate region 
to cancel the leading errors in the other two approximations in the regions where they are 
not valid. The second and third term together are dominated by values of of order A^ 
because the numerator 1 — vanishes by assumption for smaller values of p"^. 

The first term is independent of the cutoff and gives —^/B^/{8^/2Tc). The last two terms 
give D{E) a correction behaving as 1/A, with a coefficient that depends on how Ua behaves 
for p ~ A. For the Gaussian form of Ua that we have adopted, the A-dependent term in 
D{E) is 

2(27r)3/2A • ^ 
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If we want to control this cutoff dependence to renormalize the scattering matrix, we can 
vary the coefficient of 1/A by introducing an extra couphng parameter into the definition 
of t/A, which thereby would become a new coupling parameter in the potential V. There 
are many ways to introduce this parameter, all of which give higher order corrections to 
D{E) as well as a 1/A correction. We have chosen-quite arbitrarily-to use the form already 
specified for including the extra coupling constant /12(A). Universality should insure that 
this arbitrary choice is as good as any other. 

Now one finds that the A-dependent correction to D{E) takes the form 

[IQ - 8h,{A) - hl{A)]iE + B2) 

32(27r)3/2A ■ ^ 

Holding this term constant yields the renormalization group equation: 

A J^[16 - 8/i2(A) - hl{A)] = [16 - 8/i2(A) - hl{A)] . (27) 

We note that this equation has fixed points at 

h2{A) = h; = ± V2) . (28) 

At the fixed points, there is no correction to D[E) in order 1/A. The fact that ^ Q 
indicates that the operator it multiplies is not really an irrelevant operator. To isolate an 
irrelevant operator we must consider deviations of h2 from 

We close this section with a brief discussion about what happens if one uses the additional 
coupling /12 to change the effective range r^. (see Eq. (7)) from its zero continuum value to 
a finite value in order to reproduce data and use this Hamiltonian in effective field theory. 
As we have shown above, for theories with a continuum limit re = 0; but to use this as an 
effective field theory we want to be able to adjust rg. 

We can use Eq.®, with U{p) = (1 + /i2 / A^) exp(— p^/A^) to find equations for G2 and 
h2 that allow us to fix a and re at whatever values we please: 

^ 1 _ A{?,hl + Sh2 + 16) ^ 

G2 87ra 128v^7r3/2 ^^""^ 

1,9/2 2\, /2 1 4\ 

^2+H^--T^ /^2+ -r^ + -re--=- =0. (30) 



4v^A ' Vv^A aA^ \aA^ 2 ^A) 

We can fix a and re at cutoff independent experimental values by using cutoff dependent 
values for G2 and h2. In addition, the error in other A-dependent observable quantities will 
now be of 0(^2 /A^). 

Consider what happens to /i2 as A ^ 00. In this limit the equation for h2 simplifies to 

hi + 8/i2 + 2V2^Are - 16 = 0. (31) 

The solution is 



/i2 = -4 ± V32 - 2V27rAre. (32) 
This result makes it clear that we can use this model as an effective field theory only for 



cutoffs that satisfy Are < 8^2/71. Furthermore, we conclude that if we choose re 7^ it is 
not possible to completely remove cutoff dependence; because if this were possible, we would 
be able to let A ^ 00. Of course, this is not a serious limitation as long as re/a -C 1; but if 
this constraint is not satisfied by the data we wish to model, there is no reason to believe 
that we can use this simple effective field theory. 
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IV. LIMIT CYCLE FOR THE THREE-BODY COUPLING 



We now provide the initial analysis showing that a single coupling constant can approach 
either a fixed point or a limit cycle as A goes to infinity, with no other possibilities, and then 
observe that the three-body coupling of Bedaque et al. S 0] exhibits limit cycle behavior. 
We show that an infinite sequence of bound states approaching zero energy geometrically is 
a natural consequence of a limit cycle. 

Let Gz{K) be a three-body coupling constant satisfying the renormalization group equa- 
tion 

A^ = /5g3(G3). (33) 

We assume that I3g-j,{Gz) is differentiable, as it turns out to be for the example of Bedaque 
et al.. We observe that G3(A) cannot oscillate in value, if the oscillations have either a finite 
maximum or a finite minimum, because jSc-^iG^) would have to vanish at either a maximum 
or a minimum. Such vanishing would ensure that there are fixed points at the maximum or 
minimum, and G3(A) would approach one of these fixed points as its limit for A going to 
infinity, rather than oscillate. 

There are only two possibilities left for G3(A): 

1. It can change only in one direction, always increasing or always decreasing, in which 
case it has a limiting value - either a finite fixed point or ±cxd. 

2. It can oscillate but only by jumping from -|-oo to — oo or vice versa. 

The formula for G3(A) can be obtained from the -f^(A) of Bedaque et al. by multiplying 
H{h) with minus one, 

^ /AN ^ sin(goln(A/AJ - arctg(l/so)) 

' sin(soln(A/A,) + arctg(l/so)) ' ^ ' 

where A^ is a dimensionful constant that must be chosen and Sq ~ 1.00624. The additional 
minus sign is necessary to translate from the Lagrangian to the Hamiltonian formalism we 
are using. G3(A) increases while finite, until it reaches +oo; it jumps to — oo and then 
increases again. We will show in a later stage that 6*3 (A) must increase rather than decrease 
as A increases, except for the discrete jumps. 

Now we use the renormalization group equation to show that in case 2, G'3(A) must 
exhibit limit cycle behavior. We introduce a new independent variable: 

t = ln(A/Ao) , (35) 

where Aq is arbitrary, so that the RG equation will be translation invariant in t. The 
renormalization group equation now reads 

^ = /5a3(G3). (36) 

But now by assumption the solution G^{t) passes from —00 to +00 between two values of 
t, say ti and ^2- However, since the renormalization group equation is translation invariant 
in t, 6*3 (t) must repeat itself with a period t2 — ti, where ^2 > ^i- 
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Now we look at the structure of the three-body bound state spectrum. We consider the 
case that G2 is G2 so that B2 is zero. In this case there are an infinite set of three-body 
bound states with energies of the form 0, 01 : 

Bjp^ =A'F4Gs{A)]. (37) 

The factor provides the energy scale; the functions F„ provide dimensionless coefficients 
that depend on the three-body coupling parameter Gs{A). The ground state is labeled by 
n = 0. When G^IA) is near —00, Fq is near one (as will be shown later). Then as G^IA) 
increases towards +00, the i^„'s increase towards zero. As part of keeping low energy physics 
unchanged as A increases, the 's for large n do not change. However, when 6*3 (A) jumps 
from +00 to —00, a new ground state emerges and what was F„ now becomes equal to 
instead: 

F„(+oo) = F„+i(-oo) . (38) 

Another observation is this. Each time G3(A) passes through zero, the F^s must take 
on the same set of values. Thus suppose that G3(Ai) is zero. Then consider what happens 
as A increases enough so that G3(A) increases to +00 (with A now A2), jumps to —00, and 
then increases to zero again, at A equal to A3. For large n, we conclude that: 

A?F„(0) =A^F„+i(0) , (39) 

due to an intermediate matching: 

A?F„(0) = A2f„(oo) = A^F„+i(-oo) = A2f„,+i(0) . (40) 

These equations establish that the ratio of successive bound state energies, for large enough 
n, is equal to the square of the ratio A3/A1; this ratio is in turn the exponential of the period 
t2 — ti of the limit cycle: 

^(fc) =71 =^^P[2(^2-ti)] , (41) 

which is valid for all n at the two-body fixed point with a zero-energy bound state and valid 
in the limit of large n when the two-body binding energy is non-zero, as we will now discuss. 



V. THREE-BODY RENORMALIZATION GROUP: A SECOND LOOK 

In this section, we consider the renormalization group equations for the three-body sector 
in the presence of a non-zero two-body energy B2 and the two two-body couplings G2 and 
/i2- We examine how G2 and /12 affect G3(A), when both are near their fixed point values 
^2(^2) ^2- We draw on results from following sections, but only of a very general kind, 
in order to arrive at somewhat surprising conclusions. 

The first conclusion justified in the following sections is that the renormalization group 
equations for G2 and h2 apply unchanged for the three-body sector, and moreover, there 
are no further irrelevant operators which affect three-body energies to relative order 1/A. 
Further couplings need to be defined only if one is interested in corrections of relative order 
InA/A^, for example. Thus the only change we need to make from our previous discussion 
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is to generalize the renormalization group equation for G'3(A) to include dependence on G2 
and 

rIC 

A-^^(3Gs{G2,G„h2) . (42) 

If G2 and /i2 are both at their fixed points, then G3(A) exhibits a limit cycle as before, which 
we denote as g^^i-^)- The limit cycle complicates the analysis for the case that G2 and h2 
are near but not at their fixed points. We can avoid these complications by setting up a 
finite difference equation linking the values of G3(A) at a discrete set of values for A that are 
all separated by integer multiples of the limit cycle period. To be specific, we let Aq be an 
initial very large cutoff, and then define Ai, A2, etc. to be a diminishing sequence of cutoffs 
separated by the limit cycle period: 

A„+i = exp[-(t2 - h)]An . (43) 

We now study departures from the exact limit cycle for (j3(A) and exact fixed points for 
G2 and /i2 in terms of small differences from them: 

Sg2„ = G2{An)-G;{h*). (44) 
Sgsn = Gs{An) - gsciK): (45) 

5/l2„ = h2{An)-h*2. (46) 



To first order in these small quantities. 



= ^ . (47) 



where a is a numerical constant and 



= ^ . (48) 

This last equation comes from Eq. (27). 
For 5g3n we have a recursion formula: 

Sg3,,+, = Sgsn + boSg2r, + &25/i2„ + bo2Sg2jh2„ + ... . (49) 

The recursion formula includes all further possible polynomial terms in the three variables 
Sg2„, Sg^^, and Sh2„, except that no further powers of 5(?3^ by itself can appear: such terms 
would be contrary to the existence of the limit cycle when 6g2„ and 6h2„ are both zero. 

The recursion formula is easily solved by iteration, to determine ^513^ in terms of Sg^^, 
B2, Sh2o, and the ratio A^/Aq. Of special interest is that in addition to powers of A„/Ao, a 
logarithm appears too, in the form of a factor of n. The factor of n is generated from the 
602 term in the recursion formula, as can easily be verified. As a result, Sg^^ is 

- ^ (na)(5h2„)VB^ 

5g^n = Sg3o + . +■■■ ■ (50) 

This simple analysis leads us to expect that if we use a single three-body interaction, 
will display residual cutoff dependence of order S2ln(A)/A^. 
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Finally, we establish over the next few sections that the eigenvalues -B3 , for large enough 
77,, are given in terms of a universal function F: 

Bt'>=AlF{g,„,g,^,h,J (51) 

with an error of order In Aq/Aq. This result enables the use of "effective field theory" methods 
applied to the three-body case, with errors of order lnA'/(A')^ where A' is the cutoff used 
in the effective Hamiltonian. 



VI. INTEGRAL EQUATIONS FOR A CUTOFF THREE-BODY SYSTEM 

With an exact renormalization of the two-body problem, we can continue to the three- 
body problem. The two-body interaction, with couplings G2 and h2, is determined and we 
will see that a three-body interaction is required to obtain a A — >■ cxo limit for the three-body 
problem 0, Q . Once again, we are interested in first obtaining the continuum limit and 
then controlling cutoff dependence using an expansion in inverse powers of the cutoff so that 
this model can be used as an effective field theory in both the two-body and the three-body 
sectors. 

One difference between our calculations and that of Bedaque et al. is that we consistently 
cut off both the two-body and three-body interactions using Gaussian cutoffs. We believe 
that this is ultimately necessary if we want to tune the theory away from the continuum 
limit, because as we have seen the continuum limit will not typically correspond even at low 
energies to the theory we wish to model {e.g., QED or QCD). The same procedure was used 
in a recent extension of effective field theory methods to the four-body problem by Platter 
et al 0. 

A Gaussian cutoff is chosen because it enables us to use very accurate numerical methods. 
In particular, we are able to obtain exponential convergence for a uniform logarithmic scale 
[23^ in the finite difference equations we ultimately solve, as is explained in detail in Mohr's 
thesis 3|. We maintain 10-12 digits of precision over about 50-70 orders of magnitude of 



momentum when solving the bound state integral equation numerically. 
In position space, the three-body bound-state equation takes the form 

-i?3V'(ri,r2,r3) = [-V? - - 

- 92 5^(ri - r2) - g2 5^(r2 - rs) - g2 S^^rs - ri) 

+ 93 5^(ri - ra) 6^{r2 - rg)] ^(ri, ra, rs). (52) 

Here, S3 is the three-body bound-state energy, and g^ is the dimensionful coupling strength 
of a three-body contact interaction which acts only when all three particles are at the same 
point. We will exchange the dimensionful couplings g2 and g^ for the dimensionless couplings 
G2 and G3 when we regulate the interaction. The form of the three-body interaction is 
somewhat arbitrary. We have chosen the product of two delta functions because it is the 
simplest one that is non-zero only when ri = r2 = rs. Other forms could be used, but the 
results will remain unchanged if universality holds. 

In order to simplify the three-body Schrodinger equation in momentum space we can 
assume without loss of generality that Pi + P2 + Ps = 0, and we can use the symmetry of 
the wave function under exchange of coordinates, 0(pi,P25P3) = 0(Pi)P3!P2)5 etc. 
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To regulate we again replace delta functions by cutoff functions of the form 

£^A(ri - ra) = y" c/r expfip ■ (ri - r2)]f/A(p) . (53) 
The cutoff form of the two-body potential is: 

ri + r2 - {r[ + r'g 



^(ri, r2, r[, r'^) = -fi'2?7A(ri - r2)t/A(ri - ra)^ 



(54) 



In general we need to use different regulator functions, U2 and t/3, in the two-body and 
three-body interactions in order to introduce separate irrelevant operators in each sector. 
For example, this allows us to use h2 as above and introduce in the three-body sector. For 
notational simplicity we will usually drop the A subscript since these couplings are always 
cutoff dependent. 

We omit details of all the transformations to momentum space (see Ref. 0]), instead 
turning directly to the Schrodinger equation in momentum space for an energy eigenvalue 
—B^. In the center-of-mass frame the Schrodinger equation is: 

-Bsip{pi,p2) = (P? + P2 + P3)</'(Pl>P2)-^2f/2(pi+P2/2)<I>(p2) (55) 
-^2t/2(p2 + Pl/2)$(pi) - ^2t/2(pi/2 - P2/2)<I>(P3) 
+^3^3(Pl)f/3(P2)f/3(P3)$l , 

where 

*(P) = ^ / rf3gf/2(q + p/2)^(q,p) , (56) 
<^>i = ^ / c^'Pi / d^P2Us{pi)U,{p2)U3ipsMpi,p2) , (57) 

and 

93 = ■ (58) 

In the following we further restrict if to be in an S-wave of the total angular momentum. 
This means that $(p) depends only on the scalar variable p, where is p^. 

The Schrodinger equation allows v^(pi,P2) to be expressed in terms of B^, the function 
$(p) and $1. In consequence, one can derive an integral equation for $(p) which has to be 
solved simultaneously with an algebraic equation for $i. The complete set of bound state 
equations is: 



D (^-|p2 _ Jo 47r2 J -I 2p2 + 2g2 + 2pqz + B. 
g2D(-lp^-B, 



3 



(59) 

DiE)^l-vf4^,^, (60) 



92 J (2ir)S 2q'-E 
471^ Jo ^ ^J-1 2p2 + 2g2 + 2pgz + 53 ' ^ ^ 
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*i = ^ /n" [I'DMrn] - 93 ^2$i, (62) 



where p ■ q = PQ'^;. The quantities $1, -Di, and all depend on and A, although 

this has not been indicated explicitly. The quantities $1, and D all depend on B2 as 

well, and D depends implicitly also on A. 

If 6*3 is zero, only the equation for ^{p) matters, and it will have solutions only when B^ 
is at an eigenvalue. If G3 is non-zero, then one can generate a solution for any value of B^, 
namely by specifying a value for the product Gs^i- This specification turns the equation 
for into an inhomogeneous equation that should pose no constraint on B-^. Once the 
equation for has been solved, the equation for $1 can be used to determine $1 and 
hence G3 as well. We will see that if we fix G*3(A) to remove cutoff dependence for one 
eigenvalue, so that the A — >^ 00 limit can be explicitly taken, the limit for a complete tower 
of eigenvalues will be determined. 

In order to precisely solve this set of equations we introduce a method that allows us to 
disentangle small, intermediate and large momentum regions. This method is an analog of 
Efimov's disentanglement of short, intermediate and long distances in position representation 



VII. UNIFORMLY VALID EXPANSION FOR A SIMPLE FUNCTION 

The next problem we address is the determination of the limiting behavior of the solution 
of the integral equation for very large A, including subdominant terms that are smaller 
by a factor of A or A^ than the leading term. We want to first fix the bound state equation 
for the infinite cutoff limit and then develop power series expansions in inverse powers of 
the cutoff for the eigenvalue and wave function. 

We want to solve an equation involving small binding energies in comparison to the 
cutoff, because we want to start in the infinite cutoff limit. It will turn out that we can 
derive the leading order equations analytically and show that they reduce to an equation for 
the low-momentum part of the wave functions that must match at intermediate momenta a 
function that we derive analytically that depends on a single angular variable. We have a 
small three-body binding energy B^ and we will allow a small two-body binding energy B2, 
so 5/A^ is a small quantity. But the momentum in the bound-state equation ranges from 
small values all the way up to infinity, with the wave function dying exponentially above the 
cutoff. 

The strategy we use is to divide the range of p into three parts: < p < y/B^, y/B^ <^ 
p <^ A, and A < p < 00. We use separate numerical computations to handle the lowest and 
highest of these ranges, while solving the middle range analytically. 

There is a formal analysis that underlies our breakdown of the problem into the three 
parts. We will introduce the formal analysis with a far simpler problem. We consider a 
simple function of three variables: rj, p, and A. The function we examine is: 

^<"-^-^'° (p+>,)(p + A) - 
We will offer a uniformly valid approximation scheme for /, valid over the whole range of p. 
The expansion converges as a power series in ///A, but we cannot naively expand because p 



15 



ranges from order rj to infinity. We will then discuss how to compute integrals over the full 
range of p,0 < p < oo, using this scheme. 

The formula is built from three expansions for / that are valid only in limited domains of 
p. Each of the limited expansions is valid in one of the three ranges for p already mentioned. 
We assume <S A, so in every domain we can expand in powers of r]/A. In addition, when 
p ;6 77 ^ A, we can also expand in powers of p/A. We refer to this expansion as fi{ri,p,A). 
In the highest range for p, p > A, we can expand instead in r]/p; we refer to this second 
expansion as //i(r^,p, A). In the middle range ?] <^ p <C A, we can use a double expansion in 
powers of both p/A and rj/p; we call this double expansion fdiv^P, A). Our uniformly valid 
formula is 

f{v,P,A) = fi{r],p,A) + fh{v,P,A)-f47],p,A) . (65) 

The simple, intuitive way to understand how this formula allows us to produce an approx- 
imation for the full function that is valid to any given order in 7] /A in every region of 
momentum is to observe that fd will equal fi when r] <^ p and it will equal fh when p ^ A. 
Thus in these regions, where // and fh are not valid, they are exactly cancelled by fd to 
whatever order we choose to work. In the intermediate region, both /; and fh are valid, so 
fd simply cancels a double-counting in this region. 

Suppose, for example, that one wants an approximation for / that is uniformly valid 
apart from errors of third order in the ratio rj/A. Then the rule is that we first construct 
truncated versions of fi{r),p,A) and fh{i],p,A), valid to this accuracy in the corresponding 
range for p. Thus, for p ~ 77, we expand / through order (p/A)^, neglecting order (p/A)^. 
The resulting truncated expansion is: 

■^'^'^'^'^^ " (p + 7?)A ~ {p + r))A^ + {p + rj)A^ ' ^^^^ 
Similarly, we expand / through order {rj/pY to give the truncated expansion for large p: 

Now we come to the crucial rule: to determine the truncated form of the double expansion 
fd{r],p,A), one first constructs the doubly expanded versions of fi{'q,p,A) and fh{ri,p,A). 

One starts from their truncated forms just given, and in the double expansions, one keeps all 
terms that arc necessary to preserve the accuracy of the double expansion through relative 
order {rj/Ay at the opposite end of the range of p. For example, the doubly expanded version 
of fi{r),p, A), which we can denote by fidiv^P^ A), is computed by expanding fi{r],p, A). We 
assume that p is of order A, and keep terms through relative order (77/A)^: 

/«(.,, P, A) = (i^ - ^ + - - ^ + ^) + - ^ + . (68) 

We hkewise expand fh{ri,p,A) through relative order {r]/Ay assuming that p is of order 77, 
giving: 

We note that the two double expansions are identical to each other, apart from a rearrange- 
ment of terms. We beheve that this is always the case when applying the rules stated above. 
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although we offer no proof of this claim. Now the rule for truncating the complete double 
expansion fd{ri,p,A) is to set it equal to these identical double truncations: 

fdiv, P, A) = fidiv, P, A) = fhdiv, P, A) . (70) 

Given these definitions for truncated expansions to include in the formula for /, does the 
formula work? 

To show that the formula works, one has to consider the three ranges of p separately. In 
the low range: p ^ rj, where / is approximated by fi{ri,p,A), the other two terms cancel 
to the required accuracy. In the high range, p ;^ A, / is approximated by fh{r],p,A) and 
the other two terms cancel to the required accuracy. Thus the key region to discuss is the 
intermediate range 1] <^ p A. What happens in the intermediate range is that both 
fi{rj,p,A) and fh{rj,p,A) have double expansions in rj/p and p/A. But fi{r],p,A) has terms 
to all orders in 7]/p while only to second order in p/A. The function fh{rj,p,A) has the 
reverse situation - an expansion to all orders in p/A but only to second order in rj/p. When 
the two double expansions are added together, the terms that are of second order or less 
in both expansion variables are double counted. The subtraction of fd{r],p,A) removes the 
doubly counted terms. As a result the approximation for / is in error due only to terms 
which are of at least third order in both expansion parameters, and the error caused by a 
product of third order terms is 

{V/P)'{p/Af = {r^/Af . (71) 

For this simple example, the exact error can be determined, as a check on the discussion 
above. We offer two intermediate results followed by the error formula: 

fi„,A)-ffy.p.A) = ^^^^-Z^^^, (72) 

Ah,P,A)-Afep,A) = X3(?TA){(^ + ?(^) ■ P3) 

The error in the approximation for / is the difference of these two expressions, namely: 

(74) 



A^{p + r]){p + A) 



This error is smaller than f by a factor of {j]/ A)^ for all momenta, as promised. Figure ^ 
shows the three components of the expansion, fd and fh- The exact function is given by 
fi + fh — fd and at this level of approximation for these components, the difference between 
the exact function and the approximation can not be seen in such a figure. 

To complete our discussion of our simple function /, we discuss a problem that arises 
when integrating the uniformly valid approximation over the whole range < p < oo. We 
assume that the integral of the exact functions converges, so clearly if we use the entire uni- 
formly valid expansion to approximate the function the integral will converge, producing an 
expansion in powers of ?7/A of the exact integral. The problem arises only if this integration 
is performed term by term; in this case some of the terms can lead to divergent integrals. 
For one example, the two most divergent terms for large p are: 

P' _ Z. (75) 
(p + r/)A3 A3' ^'^^ 
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FIG. 1: Components of the uniformly valid expansion of a simple function. Note that f = fi + 
fh — fd- For small momenta one clearly sees fh and fa will cancel, for high momenta fi and fd will 
cancel, and in the middle one set will cancel, leaving / uniformly approximated in all regions. 

each of which give rise to quadratically divergent integrals for p near infinity. The first term 
also has subdominant linear and logarithmic divergences when it is expanded in powers of 
A/p. All these divergences are cancelled: the quadratic divergences cancel between these 
two terms, and the subdominant divergences are cancelled by other terms from f^. 

The divergences must cancel exactly, but it is convenient for analytic work to have a 
procedure that makes the integral of each term separately finite, without changing the 
integral of the sum of all terms. We can accomplish this goal by introducing a standardized 
subtraction to apply to each divergent integral to make it finite. For example, we could 
suggest the following set of standardized subtractions to make all integrals finite at their 
upper limit: 

1. Positive powers of p and constant terms independent of p are dropped completely, e.g., 

is replaced by for all n > 0. 

2. A logarithmic divergent 1/p term, integrated to infinity, is subtracted out from A to 
infinity. 

According to these rules the p/A^ term would be dropped completely, positive powers of 
p for p ^ A would be subtracted from the p^/ ((p + r))A'^) term, and the 1/p component of 
this term would be subtracted out only for p > A. To illustrate how this rule works, we 
exhibit the subtracted form for the integral over the first term: 
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Note that a rule subtracting 1/p terms for all p would not work, because they would introduce 
unwanted divergences for p ^ 0. These rules do not change the overall integration of the 
uniformly valid approximation because all the subtractions cancel out. 
A similar set of rules can be derived for handling divergences for — > 0: 

1. Negative powers of p other than 1/p are dropped completely. 

2. A logarithmically divergent 1/p term is subtracted out only for < p < 77. We use 
7] rather than A because 1/p divergences for small p typically emerge from functions 
involving 77 rather than A, and this form for the subtraction ensures that such functions 
do not acquire an unwanted dependence on A as well. 

The application of this scheme is easy when given a simple function, but its real power 
is only seen when solving an integral equation in which the function is not known. Even 
though the function is not known, we can solve for the pieces of its uniformly valid expansion 
and this is our procedure for obtaining an expansion of and the wave function in powers 
of B/A^. 

In order to solve the complete set of equations, we also need an expansion for functions 
that depend on two momenta, p and q, each of which range from to 00. We do not present 
a full analysis of this expansion here, because it requires eleven regions. A full discussion is 
found in Mohr's thesis 3 • 



VIII. LEADING ORDER THREE-BODY EQUATION FOR LARGE A 

We are now prepared to solve the three-body equations in the limit A ^ 00. To do this, 
we first need to find the leading-order approximation to the full integral equations using the 
method of uniformly valid expansions. We assume that all binding energies are much less 
than the cutoff and isolate approximations that are valid in various regions of momenta. It 
is conceptually straightforward to extend the calculation to systematically include higher- 
order corrections, but these will not affect the continuum limit that we isolate as a starting 
point. We will provide only those details that are required to quahtatively understand the 
results and once again refer the reader to Mohr's thesis for details. 

We start with a few definitions that will simplify notation somewhat: 

V2 ^ /b^, (77) 

r/3 = v^. (78) 

Using ri2 and 773 in our equations allows us to deal with quantities that have the same 
dimension as the momentum variable. It also identifies their role in the expansion as the 
parameter rj from the previous chapter. In fact, we will occasionally use rj to generically 
refer to 772 and 773 in places where either is a valid alternative. 

We use dimensionless couplings, G2 = Ag2 and G3 = g^. Further, we note that G3 
enters the bound-state equation only as a multiplicative constant that is easier to work with: 

S ^ (79) 

6 will be found to remain finite for all values of the cutoff, and at the end of the calculation 
we can readily compute $1 to isolate G3 and display the limit cycle. We will find that $1 is 
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independent of to leading order, which is consistent with the fact that we have only one 
coupling G3 with which we must renormalize the entire bound state spectrum. 

Another dimensionless quantity, and perhaps the most important, is the replacement of 
the pseudo-wavef unction with 



Like ?72 and rj^, its role in the expansion is easily seen to be the same as the function / used 
in the previous chapter. Like G2 and G3, the dimensionless nature of / will simplify future 
power counting. In addition, the p"^ term in the factor p"^ + t]'^ — rjl ensures that / tends to 
be of 0{1) throughout the entire range of p while the r/l — rjl term makes / less prone to 
large numerical fluctuations than $ when B2 — B^.^ 

With these new definitions, the three-body bound-state equation ()59|) now takes the form 



The process of expanding Eq. (jHT|) is done in a series of steps. Because the 6 term 
is simply added to the integral term, each can be expanded individually and then added 
together at the end. In addition, these terms are composed of other quantities like Di{p), 
G2, etc. which have their own expansions in terms of t]/A. The goal is to split the equation 
for / into separate equations for //, fa, and fh, isolating and solving the leading order 
equations to produce a uniformly valid approximation for the wave function. The details 
of the calculation are not very enlightening^, so we simply list the complete set of leading 
order equations that must be solved. 

Pseudo-Wavefunctions: 

We are able to precisely compute a cutoff integral over one variable of the three-body 
wave function. The fundamental equations for the calculation of these pseudo-wavefunctions 
are: 



/(^2,^3,P,A) = (p'^ + V3-V2)^ip)■ 



m 




- 6 



{y^ + vl-ril)Di{v) 
G2AW {-r]l - Ip^) 



(81) 




(82) 




(83) 



^ We cannot prove thes^_at;atements a priori and must verify them numerically. 
^ See chapter 5 of Ref. 14]. 
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/ (p A) = P\ r dq C dz ^ ^'^) (P + ^ , 

jm\F. ) 2nWho{p,A)Jo J^i 2p^ + 2q^ + 2pqz ^' 

3hl + 8/i2 + 16 /ZWPiA)\ (p^\ 



The intermediate range can be determined analytically |24l |: 



fM = Acos(soln(|)+0), (85) 



where we can choose A = 1 because the bound state equations are homogeneous, ^ is a phase 
determined by boundary conditions, and sq is the real, positive solution to the equation 



—pz smh = Sn cosh 

V 6 y V 2 . 



(86) 



The value of sq to the precision we want is 1.006237825102. 

The A dependence in Eq. (85) may be a bit misleading since Eq. (83) does not contain 
it. The cutoff is required because the argument of the logarithm is dimensionless. Making 
this choice is a matter of preference, and it could just as easily have been 773. We discuss 
this renormalization prescription below; but first we need to show how choosing two binding 
energies uniquely determines one and only one bound state spectrum. 

Clearly the choice of B2 can be made independently of any modeling in the three-body 
system, so we take this and -B3 as input for Eq. (85). We assume that > B2 so that a 
stable three-body bound state exists. This determines fio{p), which must map exactly on 
to fdoip) for P ^ Vs- Once A is chosen, 9 in f^o is determined by the phase of fio, which 
is in turn determined by both B2 and B-^. Finally, fho{p) must map onto fdo{p) for p <^ A. 
The solution requires a phase and this fixes it. Given the phase in fhoip), ^1 is determined, 
which in turn yields G3. 

This explains how one three-body bound state is determined. The remaining allowed 
values of Bs must produce the same phase 9 as is produced by the value that is fixed to 
determine 9. This produces both the remaining spectrum and the pseudo-wavefunction. 

The remaining equations we require are: 



Dlo{V2,V3,P) 

DM 

Dm{p,k) 



|p2 + ^2 _ 



8^ 



TT 



V3 



P 



167r ' 

P 



(87) 
(88) 



v^vr(4A^-3/.,j>f (l-Erf 



2567r2A4 ^ 

+ 3V27r /iaAp ((8 + /12) A^ - 3/i2P^ 



exp 



•3/~ 
2A2 



(89) 
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^ (j^. ^ A (6 (hi + 4/13 + 12) + h2 (5hl + 12/13 + 12)) 

^ 576^/3 7r3/2 ' ^ ^ 

^ ^ A (6 jhl + Ah, + 12) + (5/^1 + I2h, + 12)) 
n r 1 2 /-^ - ^2(q+b)^3(g)^3(p)^3(p + q) 



^2 = Ag2, (93) 
128^2 7r-V2 

- 3/.i + 8/.2 + 16' ^^^^ 



^3 = AV3, (95) 
^ G,oXo-D2,oVA^' ^^^^ 



3 /"^^ 

^0 = dpDiho{p,A)fho{p,A), (97) 



= 8^^ io ^ io L 2p^ + 2g^ + 2pg. • 

At this point we simply need to solve these equations numerically. We demonstrate that 
we achieve 10-12 digits of precision. 

IX. DISCRETIZATION AND PRECISION 

Because a closed-form solution for either /;o or fho is unknown, we must resort to numer- 
ical calculations of these functions as well as any energies or couplings. A common method 
for solving an integral equation involves changing the integration into a sum over discrete 
points. The integral equation then becomes a matrix equation easily solved by standard 
methods. 

While this may appear straightforward, there are several practical issues to consider. 
For example, new limits on the integral equation must be determined. It is impossible to 
numerically integrate to infinity, so a suitable upper bound must be chosen. In our case. 



22 



we use a logarithmic integration scale, which is critical in almost all problems amenable 
to renormalization, so a new lower bound to replace must also be determined. How 
the discrete points are chosen must be carefully considered. Fortunately, for our problem 
exponential convergence with the number of points can be achieved, as discussed below. 

Each choice is a compromise between accuracy and size. For example, we could choose 
an upper limit that is extremely large (ensuring that we are "close" to infinity) and discrete 
points that are closely spaced (minimizing errors in the sum). The trade-off is that using 
more points requires using a larger matrix. Using N discrete points will result in a,n N x N 
matrix with elements. If all numbers are double precision decimals, even 8000 points 
would be enough to overwhelm a computer with 512 MB of memory. This does not even 
take into account the time needed to process such a matrix. Obviously, the goal is to obtain 
the desired accuracy with a minimal number of points. We will discuss a few methods that 
drastically reduce the number of points we need. 

In the following sections, we will assume that our goal is about 12 digits of accuracy. 
This high accuracy may not be necessary for most leading order calculations, but it is 
essential when studying leading corrections. Besides directly obtaining the equations for 
the A — > oo limit, one principal reason for expanding the three-body equation in powers of 
rj/A is to analyze the cutoff dependence. If we are attempting to study this behavior, we 
must be certain that our numerical errors are not larger than the corrections being studied; 
otherwise there is no way to distinguish the small corrections from the numerical "noise." 



A. Transition to Mid-Momentum Function 

One way to limit the size of the matrix is to limit the range over which the function must 
be integrated. We know that fio and fho approach fao ioi p ^ rj and p <^ A respectively. 
This allows us to replace either function with cos (sq ln(p/A) + 9) in the appropriate range. 
The point at which we can make the switch is determined by the accuracy we desire. These 
hmits are derived for the case of 12 digits of accuracy. 

For p T), the equation for fio can be written as 

r I . 4(1 + ril/p'^ - t]IIp^) g 



X 



a/Stt (^l + (2r7i)/(3p2) - h q^ + - m 

fio{V2,V3,q)- (99) 



\^ ( P^ + q^+pq \ vl vl 

\p2 ^q2_pqj 2p2 + 2^2 ^ 2pq 2p^ + 2q^ - 2pq 



Here we have treated rj/p as a small quantity and perturbatively expanded all factors. As 
long as both r}2/p and ris/p are less than 10"^^, this equation will match the one for fao to 
12 digits. Of course, 773 must be greater than 772 since we are considering only stable bound 



(in- 



states. This means that p ~ 10 i]^ sets the limit above which fiQ can be replaced by /, 
In practice, we must have enough data points above this limit to ensure that our cosine fit 
is accurate to 12 digits also. Therefore, we will use an actual limit of p = 10^^ 773. 

Similarly, we can expand the equation for fho in the region p <^ A. Notice that in 
this region Diho{p,A) approaches Di(io(A) which is proportional to A, and the function 
Dfio{p, A) becomes equal to D^oip) = \/3p/167r. In the fho integral equation, the momentum- 
dependent part of the three-body interaction becomes 
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Since the leading order mid-momentum equation has no three-body interaction term, the 
above term will equal zero to 12 digits if we choose p ~ 10~^^ A. 
The integral part for f^Q looks like 



P 



2nWf,o{p,A) Jo 



dq 



1 ^2(q+b)^2(p + |q) 
2p2 + 2q^ + 2pqz 



-1 



;ioi) 



We have already stated that Dho(p, A) approaches Ddo{p), but more importantly, it ap- 
proaches like 

' VS {hi + 8h2 - 16) \ p_' 



Dno{p,A)'-^DM 



1 



271 



(102) 



It will therefore equal D^o to 12 digits if p ~ 10^ A. This limit also applies to the integrand 
itself, so we may replace fhQ with f^o for values of p less than 0(10^^^ A). In practice however, 
we use a limit of p = 10~^^ A to ensure that we have enough points below this region to fit 
the cosine behavior. 



B. New Integration Limits 

For the case of fho, we have limited the range of integration to be 10~^^ A to oo. (Below 
this range, we use fdo-) Naturally, we cannot integrate to infinity and instead must find a 
new hmit to replace it. Let us call this hmit A. We choose A such that 



27r'^Dho{p, A) Jx J -I 2p2 + 2q^ + 2pqz 

The exponentials in U2 suggest that the integrand should die off quickly, allowing us to 
make an initial guess for A using 

e-^VA^ ^ 10-12. (104) 

This gives an initial value of A ~ 5.25 A. However, the double integral makes the analysis 
harder since we cannot determine the exact behavior. We must resort to numerical com- 
putation, and some sample calculations reveal that a limit of A = 10 A is sufficient for our 

purposes. 

From 10^^'' A down to 0, fhQ is replaced by fdo- We would like to replace the limit with 
a larger value that still maintains our desired accuracy. Even though we know the analytic 
solution for /^o, narrowing the range of integration will reduce our computational effort. 
Call this new lower limit e, which is chosen so that 



p2 ,.c rl [/2 ( q + |p) U2 ( p + |q 



We assume that e ^ 10~^^ A < p. To within 12 digits of accuracy, the z integration can be 
written as 

"1 ^^^2(q+b)f^2(p + h) _ /I ^^t/2(p/2)f/2(p) 



J-1 ^ 2»2 + 2o2 + 2«(/^ 7-1 



2p2 + 2g2 + 2pqz J-1 2p2 -|- 2pqz 

^t'.(p/2)f/,(p) 1.(^1. (106) 
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Since q p, the logarithm can be approximated as 2q/p. Our constraint for e now becomes 



1 



27,Wf,o{p,A) 



f/2 (p/2) U2 (p) r dq fM < OilQ-^^) 
Jo 



(107) 



The values of p are of the same order as A, so we expect the value of U2 to be 0{1). 
The function fdo is also 0{1), so it is simply replaced by 1 in this approximation. This 
leaves an integral with a value of e. When combined with DhQ{p,A) ~ 0{p), we find that 
e/p < O{10~^^). The smallest value for p is 10"^'' A, implying that e ^ 10"^^ A. In practice, 
this value is sufficient for 12 digits of accuracy. 

Having replaced the limits for fho, we move on to fio- Again, we must find a finite upper 
limit to substitute for infinity. For p > 10^^ 773, fio is replaced by fdo- This new limit. A, is 
determined by the condition 



(p^ + r]l- r]l) 
A7r^pDio{r]2,V3,P) 



dq 



In 



' p2 _|_ g2 _|_ pg,^ 



^ {p^ + ril-ril) 


r°° dq / 2p 


'in'^pDio{r]2,r]3,p) 
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q \q 


{p^ + vi-vD 


POO 


% fdo{q) 


27r'^ Dio{rj2,r]3,p) -> 


'a 


(p2 + nl- fill) 


POO 


dq 


271^ Dio{ri2,ris,p) j 


'a 




{p^ + r]l- r/l) 




) < C(10 


27r2 Dio{r]2,r]3,p) 





fdo{q) 



-12 



(108) 



where we have used the fact that q ^ p rj-^. For values of p much larger than 773, the 
term in Eq. ()108p is proportional to p/\. The largest value p can obtain is 10^^ 773, implying 
A = IQp'^ 7]^. Numerical calculations verify that this limit is sufficient to assure the desired 
accuracy. 

Finally, we must replace the lower limit for /^q with a non-zero value e that we assume 
to be much smaller than 773. Our requirement is that 



(p2 + _ ^2) 



47r2p(r^| -i]^)Dio{r]2,ri3,p) Jo 



dqq In 



' rjl + 2p^ + 2pg' 
ril + 2p^-2pq 



fio{V2,V3,q) <O{10 



-12\ 



(109) 



Expanding the logarithm to 0{q) yields 

(p2 + _ ^2) 



T^'^ivi - vDijll + V)Ao(^2,^3,P) -^0 



dqq^ fioiV2,V3,q) < 0(10 



-12\ 



;iio) 



For small values of q, we will find that fiQ is approximately constant and of 0{1). Therefore, 
the integral is roughly equal to e^. If p <^ 773 or p ~ rj^, then Eq. ()110|) is 0{e^/ri^). If 
p ^ f]3, then it is 0{e^/pr]^) < 0{e' /rf). This seems to imply that a value of e ~ 10"^ 773 is 
adequate. 

Unfortunately, using this value of e will result in poor accuracy when 772 — 773. This is 
a result of the [rjl — rjl) term in the denominator of Eq. ()110|) . Originally, this was our 
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approximation to the term (g^ + vi ~ V2) integral equation. When the energies are 

nearly equal, our approximation needs to be g^. The condition on e should now become 



2f 2 !o^2?n ( V rdqUv2,V3,q) < 0(10-''). (Ill) 

The integral is roughly equal to e, and the entire term is 0{e/ri). This means we must use 
the lower limit e = 10~^^?73. 

These derivations are general enough to determine the appropriate limits for other cases 
of desired accuracy. Keep in mind however that the limits must always be tested numerically 
to ensure that they are indeed sufficient. 



C. Discrete Point Spacing 

Now that we have limits in place, we must choose the discrete values of q within these 
limits at which to evaluate our functions. We employ discrete points that are equally spaced 
on a logarithmic scale. There are two main reasons for making this choice. 

First, we have already seen that fdo = cos (sq ln(p/A) + 9) is periodic on a logarithmic 
scale. In fact, other functions have similar behavior, including G3. It makes sense that a 
logarithmic spacing is suited to capturing the behavior of the system. 

Second, and far more important, by choosing points in this manner an integration of q 
from to 00 becomes an integration of ln(g) from — cxo to 00. By equally spacing points 
on this log scale, we can achieve convergence that improves exponentially with the spacing 
[isj ]. as opposed to the more typical power law convergence. This drastically reduces the 
number of points needed to achieve our desired accuracy. For instance, to cover the range 
of fio from 10~^^773 to 10^^ 173, and maintain 12 digits of accuracy, we can space our points 
by Pn+i = Pn^^'"^- This requires only about 300 points. Exponential convergence with 
the number of point is extremely powerful and it is not widely appreciated that this is ever 
possible. A simple demonstration (with minor errors) can be found in an appendix in Mohr's 
thesis fl4^. Figure |21 illustrates this exponential convergence which eventually terminates due 
to roundoff errors. 



X. DISCRETIZED EQUATIONS 

The discretized integral equations are shown below. Here, A represents the momentum 
spacing, and the momentum values are related by the equation Pn+i = Pn^^- The identity 
matrix is represented by Inmi and the ceiling function represented by \x\ returns an integer 
value n such that {n — 1) < x < n. 

Low-Momentum: These are easily derived from Eq. (83). The value of n for p„ ranges 
from to Nmid, while the value of m for ranges from to Nmax- These are defined as: 
po = 10-12 r/3, N^,d = rin(1027)/Al, and K^ax = rin(1039)/Al . 

{Mnm - Inrn) flo{Pm) = K, (112) 

m=0 

^ ^ ^ qI (pI + Vi- Vl) ( Vi + + 2ql + 2pnqm \ 

47r2Ao(r?2, Pn) Pn{qL + ri- vD \vl + '^pI + - 2pnqm) ' 
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FIG. 2: Relative error in the phase appearing in pseudo-wave functions as a function of the number 
of points used in the discretized approximation for several pairs of two-body and three-body binding 
energies. 



br,. = 



9m 



X COS ( So In ( — I + ^ 



(114) 



High-Momentum: The value of n for p„ ranges from to Nmax-, while the value of m for 
ranges from -N^^^ to N^a.x- These are defined as: po = 10"^'^ A, N^aax = Rn(10^^)/A] , 
and Nmin = rin(10i2)/A] . 



X] (-^nm - 4m) fho{Pm) = On + (^0 &n, 



m=0 
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X K{pn, Qm) COS (^Sq In (^^^ + 9^ , (117) 
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(119) 
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The hyperbolic tangent discretization we use for the angular integral, with 401 angular 
points, has been numerically verified to exceed our needs over the entire range of binding 
energies we explore. There may be equally good or better choices. 

XI. ANALYTIC AND NUMERICAL RESULTS 

A computer program to solve the discretized integral equations can be implemented with 
the help of some basic numerical algorithms. Once completed, it can be used to generate 
highly accurate results and analyze the three-body system. Nonetheless, attempts to extract 
analytic results should not be overlooked. Even some of the most general properties of the 
integral equations allow us to draw conclusions about the behavior of the system. 

We begin this chapter with analytic results obtained from studying the leading order 
integral equations. These results include statements about the cutoff dependence of bound- 
state energies and the phase for f^Q- A proof for the cyclic behavior of 5q is given, which is 
then used to infer similar behavior for G^. 

Following the analytic results is a section containing numerical solutions to the integral 
equations. Here we examine behavior that cannot be determined analytically. Solutions for 
the functions //o and fho are shown, and the cutoff dependence of G3 is calculated. Some 
relations that are proven analytically are also verified numerically. 

In Section IVIIH we saw that Eq. ()82j) contains no A dependence, making the function 
//o(p) cutoff- independent. However, for values of p ^ 773, the solution must behave like 
fdo{p) = COS (so In (p/A) + 9), which explicitly has the cutoff in it. The only way these two 
equations can be reconciled is if 6 contains cutoff dependence of the form sq ln(A). Any 
remaining part of 6 must be a function of 772 and rj^. Therefore, we can write 



where ^ is a dimensionless function of the ratio ?72/'73- This relation holds for any values of 
ri2 and rj^, including all 773 values corresponding to multiple bound states with the same 772- 
Using 773 in the ratio with A is a matter of choice. Any quantity composed of 772 and 773 with 
the same dimension as A would work just as well; however, we need to allow 772 = 0, so 772 
alone is a poor choice. 

While Eq. fjl2H) is quite simple, it has many interesting consequences. As we mentioned 
earlier, two different three-body bound states in the same spectrum must have the same 
phase 9. Suppose that we choose some fixed values for 772 and 773 and make them cutoff- 
independent by choosing the appropriate A dependence for G2 and G3. This might be 
desirable if we are trying to match those energies to experimental data. The phase for the 
solution in this case is 



^ = So In (A/7/3) + 9 (7/2/773) 



(121) 



9a = So In (A/773) + 9 (772/773) , 



(122) 
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for some givpn cutoff A. If the cutoff is changed, G2, G3, and 9a will all change with it, but 
772, Vs, and ^(772/^3) will not. 

Imagine now that we find a second three-body bound-state solution with the same phase. 
Let us call its energy 773. The phase for this solution is 

eA^soin{A/fj3) + e{m/f)3), (123) 

which must be equal to by assumption. This results in the relation 

So In (773/773) = e (772/773) - e (772/773) . (124) 

If the cutoff is now changed to a new value A', the original data gives a phase of 

9a> = So In (A7773) + 9 (772/773) . (125) 

The question is whether 773 is still a valid solution. Its new phase is 

9a' = So In (A7773) + ^(772/773) 

= So In (A7773) + [so In (773/773) + 9 (772/773) 

= So In (A7773) + ^ (772/773) 

= 9a'. (126) 

Since the phases still match, 773 is still a bound-state solution. This remains true for any cut- 
off, implying that 773 is also cutoff-independent like 773. Of course, the same statement applies 
to any other three-body bound state making the entire spectrum completely independent of 
A. Such behavior should come as no surprise since the leading order equations represent the 
A ^ 00 limit. Keep in mind that this is true for any other physical quantity, but does not 
apply to the couplings. Obviously, it is the cutoff dependence of the couplings that enables 
the bound states to be cutoff independent. We have shown that a single three-body contact 
interaction allows us to renormalize the entire three-body bound-state spectrum. 

Just as fio has no dependence on A, neither does it have any dependence on h2- As 
a consequence, /12 has no effect at leading order on the binding energies or other physical 
quantities. It does appear in the equation for fn however, showing that it is needed when 
considering first order corrections. Because we work only to leading order, we use /12 = 
below unless stated otherwise. 

We now turn to the equation for fho{p)' 



p2 fl 



U2 (q + h) U2 fp + k) 

_ ?,hl + 8/i2 + 16 ( DimjpA) 
' 128^2 7r3/2 \dUp,^) 




(127) 



The dependence on p and A has been explicitly shown. Since fho is dimensionless, only the 
ratio p/ K can occur in the function. Furthermore, there is a dependence upon h2 and 9. 
The /i2 dependence comes from its appearance in the function U2{p), either directly in the 
integral or indirectly in Dho and Diho- The 9 dependence is a result of fho approaching fao 
for p < A. 
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Assuming that A is held fixed, choosing a value for So will uniquely determine the phase. 
Thus, we can view the phase as a function of the coupling, 6 {60). Conversely, choosing a 
phase determines the coupling, so it is equally valid to treat the coupling as a function of 
the phase, 50(6*). The coupling 60 can also have a dependence upon h2 but not upon A. 
The reason is that 60 is dimensionless and there is no other quantity available to form a 
dimensionless ratio with A. 

With this in mind, consider a solution to Eq. (jl27|) with a coupling of So{6) and a function 
fho{p) that behaves asymptotically as cos (sq In (p/A) + 9) for small p. For a phase of 9 + 27i, 
the cosine behavior will remain unchanged. We may therefore conclude that 

6o{e + 27r) = 6o{e), (128) 

showing that the coupling exhibits periodic behavior. If we apply the operator —d'^/dO'^ to 
both sides of p27|) . we find once again that the asymptotic cosine behavior is the same. 
This implies 

- ^ = m- (129) 

The conclusion to be drawn is that the coupling may be written as a cosine function with 
some amplitude A and phase 0. These two parameters will contain any /i2 dependence that 
9 may possess, so we shall write 

So{9) = A{h2) cos {9 + 4>{h2)) . (130) 

The equation relating 60 and G3, 

^3,0 = 7^ — -p^ — TTaI' (^^^) 

G2,0-io - 1^2,0 cio/ A* 

shows that whenever Sq is zero, so is G^,.^ Since two adjacent zeros of Sq occur when the 
cutoffs satisfy 

9a' -9a = so In (^j^ = IT, (132) 

the adjacent zeros of G3(A) should be spaced by a cutoff factor of A'/A = exp(7r/so) — 
22.69438259536. This suggests that G3(A) may also possess cyclic behavior, but this must 
be verified numerically, which is done below. 

We begin our numerical investigation by considering solutions for the functions fio and 
fho which lead to an approximation for the complete function f{p). Next, the cutoff in- 
dependence of the three-body spectrum is verified, followed by an analysis^f the coupling 
constants. These numerical results are confirmed by Wilson's calculations 



A. Pseudo-Wave Functions 

Figure 01 shows the numerical solution for fio{p) in the case of B2 = 0.1, B3 = 1.0, and 
A = 10^. Notice that it is constant for small values of momentum and then takes on the 



^ The only problem that might arise is if Xq = in such a way that the ratio is non-zero, but we have found 
no parameters for which this is true. 
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f^Q(x) = cos(x + 5.684386276089572) 
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FIG. 3: Numerical solution of fio{p) shown matching on to fdo{p) for the case of B2 
and A = 10*^. 



0.1, B3 = 1.0, 



cosine behavior of f^o as p becomes large. Since the behavior of the solution to the low- 
momentum equation is determined by the ratio -B^Z-B-^, this figure is representative of all 
solutions, as demonstrated explicitly in Mohr's thesis 3 • 

This binding energy leads to 6 = 5.684386276089572..., and with A = 10® we find the 
numerical solution for the function fho{p) shown in Figure 0] This is the high-momentum 
function corresponding to the low-momentum function in Figure El For all fho solutions, we 
see the cosine behavior for p <^ A and a suppression of large momentum values when p > A. 
This suppression is an effect of the Gaussian behavior of 1/2- 

Combining fiQ, f^o and fho, we see the overall leading order behavior for the function 
f{p). Figure El shows fio + fho — fdo for B2 = 0.1, = 1.0, and A = 10®. All solutions 
show this same qualitative behavior, with a flat plateau at low momenta going to cosine 
behavior for medium momenta and finally decaying rapidly at large momenta above the 
cutoff. Increasing the cutoff while changing the coupling so that the phase remains constant 
simply extends the cosine region, pushing the eventual decay to higher momenta. 

Several bound states typically exist for the same values of B2 and A. In the case of 
B2 = 1.0 and B3 = 1.0 for A = 10®, there exist states of energy B3 = 6.7502901502599 and 
i?3 = 1406.130393204. The solutions of fio for these energies are shown in Figure IHl AH of 
the functions are very similar. The only real difference is in the length of the initial plateau. 
As the bound-state energies become larger, the flat region becomes longer, joining the cosine 
near later peaks. For a fixed cutoff there are a finite number of solutions because the cosine 
terminates near the cutoff, but as the cutoff is extended and further periods of the cosine 
result, additional bound states appear. 
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X = Sq ln(p/A) 



FIG. 4: Numerical solution of fho{p) for the case 6 = 5.684386276089572 and A = 10^. The dashed 
line is the best-fit cosine curve that matches the low-momentum behavior. 
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FIG. 5: Numerical solution of /o(p) for the case of B2 = 0.1, S3 = 1.0, and A = 10^. 




FIG. 6: Numerical solution of fio{p) for three bound states using B2 = 1.0 and A = 10®. The 
three-body bound state energies are 1.0, 6.7502901502599 and 1406.1303932044. 

B. Three-Body Binding Energies 

Above we proved that the three-body bound-state spectrum is cutoff independent to 
leading order. We choose B2 = (for which there is a limit cycle that persists to arbitrarily 
small cutoffs) and let change with A so that the state i?3 = 1.0 is held constant. For this 
special case there are an infinite number of low- lying bound states as shown by Efimov 0, |^ . 
Two other states, one shallower and one deeper, are calculated as the cutoff changes. Since 
very small fluctuations are impossible to see in a plot. Table HI shows the calculated energies 
for several values of A. Using Efimov's result that the ratio of adjacent binding energies is 
exp(27r/so) when -B2 = 0,1^, the relative error for each calculation can be determined and 
is also given in the table. This illustrates that each energy is cutoff-independent to about 
12 digits and also matches the true value to the same accuracy. As an additional example. 
Table |n] shows the case of -B2 = -B3 = 1-0 and considers the next two deeper states as A 
changes. The binding energies are approximately 6.75029 and 1406.13. These energies have 
been previously calculated by Bedaque et al. and by Braaten, Hammer and Kusunoki 
(20I I using a method that gives at most two digits of numerical precision. Their results are 6.8 
and 1.4 x 10^, which match the results given here to within their relative errors of (9(10^^). 



C. Couplings 

Equation (jl3Up shows that the coupling 5q should have a cosine dependence on the phase 
9, which is defined by /^o = cos(so ln(p/A) + 9). Numerical data for 5q as a function of the 
phase show this behavior precisely jiU], independent of B2 and -B3. 

We suggested above that this periodic behavior should carry over to G^. Figure d displays 
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A 


^3 (Shallow) 


Error 


S3 (Deep) 


Error 


lUUUUU.UUUUUUUU 
738905.60989306 
5459815.0033144 
109663315.84284 
3631550267.4246 


0.0019416156131338 
0.0019416156131358 
0.0019416156131435 
0.0019416156131358 
0.0019416156131435 


6.7e-13 
3.2e-13 
4.3e-12 
3.2e-13 
4.3e-12 


515.03500138461 
515.03500138403 
515.03500138287 
515.03500138520 
515.03500138520 


5.2e-13 
1.6e-12 
3.9e-12 
6.0e-13 
6.0e-13 



TABLE I: Binding energies and relative errors for the next shallowest and next deepest 3-body 
bound states for B2 = 0.0, Bj, = 1.0 and various cutoffs. For this case we know the spectrum 
analytically and the determination of errors is straightforward. 



A 


B3 #1 


B3 #2 


100000.00000000 
738905.60989306 
2008553.6923187 
14841315.910257 
298095798.70417 
5987414171.5197 
44241339200.892 


6.750290150257678 
6.750290150257678 
6.750290150255419 
6.750290150268966 
6.750290150257678 
6.750290150259935 
6.750290150257678 


1406.13039320296 
1406.13039320593 
1406.13039320345 
1406.13039320345 
1406.13039320593 
1406.13039320345 
1406.13039320296 



TABLE IL Binding energies of the two next deeper 3-body bound states for B2 = = 1.0 and 
various cutoffs. 



a plot of (j3 as a function of A. This data exhibits the hmit-cycle behavior of the three-body 
coupling. As the cutoff increases, G3 becomes larger and larger, eventually diverging to 
+00. It then jumps to —00 and increases again. This limit-cycle behavior is not dependent 
upon any specific bound-state values, but the positioning of the cycle is dependent upon the 
energies. This cyclic behavior has been previously observed |y| using a sharp cutoff that 
simply discards all momenta above A and a different method for including the two-body 
interaction. 

It is shown above that h2 has no effect on the binding energies to leading order, but 
it does affect G3. In fact, Eq. (|13U|) explicitly exhibits such dependence. In Figure |H1 we 
illustrate the effect on of changing /12. The limit cycle behavior persists and the curves 
are simply distorted by the presence of non-zero h2- Again, this behavior has been studied 
for many cases 3 . 



XII. CONCLUSIONS 



We reduce the equation for S-wave bound states of three bosons interacting via attrac- 
tive short range two-body interactions to a set of coupled integral equations and develop 
numerical methods that enable us to solve for both eigenvalues and pseudo-wave functions 
with a precision of 11-12 digits. This problem must be regulated (a Gaussian cutoff here to 
achieve high numerical precision) and renormalization produces a well-defined infinite cutoff 
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FIG. 7: The dimensionless three-body couphng G3 as a function of the cutoff A for B2 = 0.5 and 
= 1.0. The limit-cycle behavior of G3 is evident from the fact that it is periodic in ln(A). 
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FIG. 8: The dimensionless three-body coupling G3 as a function of the cutoff A for B2 = 0.5 and 
B3 = 1.0 with two different values of /t2- 
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limit {i.e., a continuum limit) if a three-body contact term is introduced with its coupling 
precisely tuned to a limit cycle. 

The method of uniformly valid expansion is developed to separate widely different regions 
of momenta as a first step towards isolating and then controlling the high momentum region. 
We believe that this method will have applications in other renormalization problems since 
it addresses a generic need, but our focus is the renormalization of the three-body problem. 

Unusually high precision is required because this renormalization problem is intrinsically 
nonperturbative and must be solved numerically at this time. We must isolate effects that 
scale as mixed powers of ln(A) and and this places stringent requirements on numerical 
renormalization if one needs to accurately resolve even a few levels of detail. Irrelevant 
operators should enable us to control these sub-leading corrections with a finite number of 
couplings, if we can identify appropriate irrelevant operators. It is simplest to assume that 
we can employ irrelevant operators from free field theory {i.e., powers of derivatives acting 
on regulated delta-functions for the few-body problem) but this needs to be demonstrated 
numerically. 

Renormalization replaces an infinite geometric tower of bound states with effective in- 
teractions as the cutoff decreases. Irrelevant operators must allow us to maintain the least 
bound states of this tower as the cutoff descends. Our leading-order calculations show how 
high-momentum effects are funneled through an intermediate-momentum shell to produce 
a cutoff-independent tower of bound states that is entirely determined (with the two-body 
interaction fixed) by a single three-body coupling. The uniformly valid expansion provides a 
tool that should expose this structure at arbitrarily high levels of precision while providing 
direct insights into coupling between widely separate scales. 

The three-body bound state problem we study has a rich history and many interesting 
current applications. The continuum limit of this problem is solved with high precision and 
the groundwork is laid for the next step of modulating residual cutoff dependence so that 
this model can be applied further. 
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